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Abstract 



A neutrino mass matrix model with a bilinear form M v = k v {MoM^ M^) 2 is proposed 
within the framework of the so-called yukawaon model, which has been proposed for the 
purpose of a unified description of the lepton mixing matrix Upmns and the quark mixing 
matrix Vckm- The model has only two adjustable parameters for the PMNS mixing and 
neutrino mass ratios. (Other parameters are fixed from the observed quark and charged 
lepton mass ratios and the CKM mixing.) The model gives reasonable values sin 2 20 i2 — 0.85 
and sin 26> 23 ~ 1 and sin 2 20 13 ~ 0.09 together with R v = Am^/Am^ ~ 0.03. Our 
^v^j , prediction of the effective neutrino mass (m) in the neutrinoless double beta decay takes a 



sizable value (m) ~ 0.0034 eV. 

PCAC numbers: 11.30.Hv, 12.15.Ff, 14.60.Pq, 12.60.-i, 



1 Introduction 

Most particle physicists think that the mass spectra and mixing patterns of quarks and 
leptons, the Cabibbo-Kobayashi-Maskawa mixing matrix Vckm [I] and the Pontecorvo-Maki- 
^ . Nakagawa-Sakata mixing matrix Upmns [2], will be understood by a unified description. As 

one of such models, "yukawaon" model [31 HI [5j [6] has been proposed. The model, which is a 
kind of "flavon" model [7], has the following characteristics: (i) Yukawa coupling constants Yf 
(/ = u, d, e, • • • ) in the standard model are understood as vacuum expectation values (VEVs) 
of scalars ("yukawaon") with 3x3 components, i.e. by y/(Y/)/A, where A is an energy scale of 
the effective theory, (ii) The hierarchical structures of the effective Yukawa coupling constants 
originate only in a fundamental VEV matrix ($o}> whose hierarchical structure is ad hoc assumed 
and whose VEVs are fixed by the observed charged lepton masses. (Of course, the goal in 
the yukawaon model is to understand the VEVs of the fundamental scalar 3>o itself on the 
basis of a dynamical consideration.) In the yukawaon model, in principle, there are no family- 
number-dependent parameters except for (<3?o)- (Regrettably, in order to obtain reasonable 
values of quark mixing matrix Vckm, at present, we need a phase matrix P u (or P^) with two 
phase parameters [6]. The final goal of our model is also to remove such family dependent 
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parameters.) (iii) Relations among VEV matrices are obtained from SUSY vacuum conditions 
for a given superpotential under family symmetries and R charges. (It is a characteristic in the 
yukawaon model that adjustable parameters are quite few because the observed charged lepton 
mass values are used as the input values.) 

In a series of previous yukawaon models [H [6], quark and lepton mass matrices were 
given as follows: 

M e = k e $ {l + a e X 3 )$ , 
M u = M D M^ x Ml, 
M D = M e , 

M R = k R (M u M e + M e M u ) + i v term, (1.1) 

P U M U P U = k u M u M u , 

M u = A4$ (l + a u X 3 )$ , 

M d = k d [$ (1 + a d X 3 )$ + mgl] , 

where M u , M d , M u , and M e are, respectively, mass matrices of up, down quarks, light neutrinos, 
and charged leptons. Md and Mr are the Dirac and right handed Majorana neutrino mass 
matrices, respectively. Here, we introduced a common mass matrix form 

M f = kf$o(l + a f X 3 )$ . (1.2) 

(For convenience, we have dropped the notations "(" and ")" on the VEV matrix (^o))- X 3 , 1 
and P u are also VEV matrices of other scalar fields. The matrices $o> 1 an d P u are defined 
by 
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and P u = diag(e _l< ^ 1 , e - ^ 2 , 1), respectively. The form (1 + a e X 3 ) is due to a family symmetry 
breaking U(3)— > S3 [6] as we discuss later. The coefficients a/ play an essential role in obtaining 
the mass ratios and mixings, while the family- number independent coefficients kf and k' u do 
not. The values of (xi,X2,x 3 ) with x\ + x\ + x\ = 1 are fixed by the observed charged lepton 
mass values under the given value of a e . In order to get reasonable fits for lepton mixing 
angles sin 2 2#23 and sin 2 2#i2 keeping reasonable fits for up-quark mass ratios, an additional 
term, "£„-term", was added in the mass matrix Mr. In this model, only the up-quark mass 
matrix M u is given by a bilinear form of P U M U P U = k u M u M u . Such a form was suggested by a 
phenomenological fact Mu ag ~ (M d ha9 ) 2 . 

In this paper, we shall present the new mass matrix model. The different points of the 
new model from the previous ones are as follows: (i) there is no ^-term and (ii) Considering 
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the correspondence between quark and lepton mass matrices (M e O M^; M v <-> M u ), we try to 
consider a form 

M y = k v M v M Vl (1.4) 

corresponding to the up-quark mass matrix form M u = k u M u M u because we assumed the form 
M d ~ M e = /c e $ (l + a e A 3 )<I>o. Here is given by 

Af„ = MdM^mI, (1.5) 

and we take 

M D = fcuSjftl + a D A 2 )$ , (1.6) 

M fl = fc K (MJ\4 + M e M„). (1.7) 

Note that the form of Mp is given by (1.6) against other mass matrix forms (1.2). The matrix 
form Xi [8] is defined by 



X 2 = I 
2 



/ 1 1 

1 1 | - I 1 ^ 

V o o o 



By this revision (1.4) with rather simple structures of yukawaon VEV structures, we can 
obtain reasonable prediction of sin 2 2#i3 ~ 0.09 together with other PMNS and CKM parameters 
and quark and lepton mass ratios. (The big drawback in the previous yukawaon model (1.1) 
was that the model could not give the observed large value [9] sin 2 2#i3 ~ 0.09.) 

In the next section, we give a new model for quark and lepton mixings and their mass ratios 
on the basis of a revised yukawaon model. In Sec. 3, we discuss parameter fitting of observed 
values only for the PMNS mixing and neutrino mass ratios because we revised the model only in 
the neutrino sector. The parameter values in the down-quark sector are effectively unchanged, 
so that we can obtain the same predictions for the down-quark mass ratios and CKM matrix 
parameters without changing the successful results in the previous paper [8] . 

2 VEV matrix relations 

We assume that a would-be Yukawa interaction is given as follows: 

W Y = ^elYVtjH, + ^H^^Hu) + V j-d ci Y^H d + ^u cl Y^H u , (2.1) 

where i = {yi,,e£) and q = (ul,cLl) are SU(2)^ doublets. Assignments of these fields to family 
symmetries U(3)xU(3)' are given in Table 1. Note that in Eq.(2.1) there are no SU(2)x singlet 
neutrinos. We have straightforwardly defined the neutrino mass matrix M v by the second term 
in Eq.(2.1). Although we denoted in Eq.(1.5) as if the matrix Mjj is a Dirac neutrino mass 
matrix, the matrix does not have a meaning of the Dirac mass matrix [see Eq.(2.3) later]. 
Under the definition of Yi (Y q ) in Eq.(2.1), the quark mixing matrix Vckm and the lepton 



3 



mixing mixing matrix Upmns are given by Vqkm = UuUd and Upmns = UjU u , respectively, 
where U f are denned by XJ]m]m s U s = Df (D f are diagonal). (Here and hereafter, sometimes, 
we denote Yi and Y q as Yf for simplify. In order to distinguish each yukawaon from others, we 
assume that Yf have different R charges from each other together with R charge conservation. 
(Of course, the R charge conservation is broken at the energy scale A.) 



Table 1: Assignments of SU(2) L xSU(3) c xU(3)xU(3)' and R charges 
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We obtain VEV matrix relations from the superpotential which is invariant under the family 
symmetries U(3)xU(3)' and R charge conservation. In the yukawaon model, the VEV matrix 
relations are phenomenological ones, and they are dependent on the R charge assignments. Since 
derivations of the VEV matrix relations are essentially similar to those in the previous papers [31 
E1E1 El [8] although the assignments of U(3) xU(3)'and R charges are different. Besides, in order 
to derive the desirable mass matrix relations, we must to consider a complicated superpotential 
form. The purpose of the present paper is not to derive those mass matrix relations uniquely, 
but it is to investigate the neutrino mass matrix M u = k u M u M u from the phenomenological 
point of view. Therefore, the derivation of the mass matrix relations are discussed in Appendix. 



In the present section, we present only the results: 

(Y e ) = M$ ) i 1 + a e XX T ) (2.2 

(X v ) = k v {Y")(Y v ), (2-3: 

(I-) = k' u (Y D )(Y R )(Y D ) = C(F d )((Vk))- 1 (^d), (2.4 

(Yd) = k D ($o) (1 + a D X T X) ($ ), (2.5 

{Y R } = k R [{Y e )(Y u ) + (Y u )(Y e }), (2.6; 

(Y u ) =k u (Y u ){Y u ), (2.7 

(Y u ) = k' u (^ )(l + a u XX T )(^), (2.8 

(P d )(Y d )(P d ) = k d ($ ) (1 + a d XX T ) + e d l, (2.9; 
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where X has phenomenologically been introduced in the previous model [8] and it has the form 



—in 

v x 




(2.10) 



The form (2.10) leads to 



((x)(x T )) = ^(x 3 u, {(x T )(x))^ = hx 2 ) lv 



(2.11) 



together with (X){X) = {X), where X 3 and X 2 is defined by Eqs. (1.3) and (1.8), respectively, 
and, for simplicity, we have put vx = 1 because we are interested only in the relative ratios 
among the family components. At present, there is no idea for the origin of this form (2.10). 
We may speculate that this form is related to a breaking pattern of U(3)xU(3)' (for example, 
discrete symmetries U(3)xU(3)' — ^xSa). In the present paper, the form (2.11) is only ad hoc 
assumption. However, as seen later, we can obtain a good fitting for the neutrino mixing angle 
sin 2 20i3 due to this assumption. 

3 Parameter fitting 

We again summarize our mass matrix model as follows: 



Y e = A; e $o(l + a e X 3 )^, 

% = KYdY^Yd, 
Y D = k D ^(l + a D e iaD X 2 )$ , 
Y R = k R (%Y U + Y u Y e ) , 



(3.1) 

(3.2) 
(3.3) 
(3.4) 
(3.5) 



Y u = k u Y u Y u , 
Y u = k' u $ {l + a u e ia -X 3 ) S% , 



P d Y a P d = fc d [$ (l + a d X 3 )*o + £>1] 



(3.6) 
(3.7) 
(3.8) 



where, for convenience, we have dropped the notations "(" and ")". In numerical calcu- 
lations, we use dimensionless expressions = diag(xi, x 2 , x 3 ) (with x\ + x\ + x\ = 1), 
Pd = diag(e _l< ^ 1 , e~ 1 ^ 2 , 1), and E = diag(l,l,l). The parameters are re-refined by Eqs. (3.1)- 
(3.8). In Eqs. (3.4) and (3.7), since we assume that the parameters a e and a d are real, while a u 
and an are complex in our Mb ^ M u and M e <-> correspondence scheme, we have denoted 
a u and as a u e tau and ooe'" 13 , respectively. 
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Table 2: Process for fitting parameters. Of course, since these parameters listed in each step 
can slightly affect predicted values listed in the other steps, we need fine tuning after the step 
5th. New parameter fitting in the present paper is started from the 5th step. 
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In this model, we have 2 parameters (a>£),aD) for neutrino sector, 4 parameters ac, £q 
and (0i,0 2 ) for down-quark mass ratios and Vqkm, and 3 parameters a e , (a u ,a u ) for charged 
lepton mass ratios and up-quark mass ratios as shown in Table 2. Especially, it is worthwhile 
noticing that the neutrino mass ratios and Upmns are described only two parameters after a e 
and (a u ,a u ) have been fixed from the observed CKM mixing and up-quark mass ratios. There 
is effectively no change in the mass matrix structures except for Y v from the previous paper [8], 
so that we can use, even in the present model too, the same parameter values for a e and (a u , a u ) 
taken in the previous study, which are given by 

a e = 7.5, {a u , a u ) = (-1.35, 7.6°), (3.9) 

Therefore, as far as PMNS mixing and neutrino mass ratios are concerned, we have only two 
free parameters {an, old) in the present neutrino mass matrix model. 
At present, the observed values [TO] are as follows: 

sin 2 2<9^ s = 0.857 ± 0.024, sin 2 20^ > 0.95, sin 2 2<9^ s = 0.098 ± 0.013, (3.10) 



RObs = (A^i)* _ (7-50 ± 0.20) xlO- 5 eV 2 _ +0 . 14 2 
Kv ~(Aml 2 )obs- ( 2 .32^ 2 ) x 10-3 eV 2 ~ t^-auJ * ™ ■ 



Since the parameters (od,0!d) are sensitive to the observables sin 2 26^ and R^ bs , we use 
the observed values of sin 2$i2 and Ry in order to fix our parameter values (o,£),oi£)). In Fig.l, 
we illustrate an allowed parameter region of {cl£),oid) obtained from the observed values [TO] 
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sin 2 26f 2 s = 0.875 ± 0.024 and R°J> S x 10 = 0.323±g;g^f. As seen in Fig.l, the observed values 
uniquely fix the parameter values {an , od) as 

(a D ,a D ) = (8.7,12°). (3.12) 

It is worthwhile noticing that the parameter values (3.12) uniquely give a prediction of sin 2 20i3 ~ 
0.09. 




Figure 1: Allowed parameter region in (od,qd) plane. The solid and dashed curves indicate 
the border and center curves of the allowed region which are obtained from the observe values 
sin 2 20^ = 0.875 ± 0.024 and R° bs x 10 = 0.324±°;°^, respectively. The dot-dashed curves 
represent contour curves of sin 2 2#i3 = 0.08, 0.09, and 0.10. 

For reference, in Fig. 2, we illustrate behaviors of sin 2 20i2 and R v versus old for the case of 
dp = 8.7. The choice = 12° gives excellent fittings to the observed values of sin 2 20i2 and 
R y simultaneously: 

sin 2 26»i2 = 0.8544, R v = 0.0331. (3.13) 
Then, we obtain our predictions for sin 2 2023 an d sin 2 2#i3 using (3.12) as follows: 

sin 2 2023 = 0.9962, sin 2 20 i3 = 0.0907, (3.14) 

which are in excellent agreement with the observed values in given in Eq.(3.10). 

The fixing of the parameters (od, od), Eq.(3.12), leads to the prediction of the CP violating 
phase parameter in the lepton sector too: 

5 e CP = 127° (J £ = 2.74 x 10" 2 ). (3.15) 

We can also predict neutrino masses: 

m u x = 0.00061 eV, m u2 = 0.00899 eV, m v3 = 0.05011 eV, (3.16) 
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Figure 2: Lepton mixing parameters sin 2 29\2, sin 2 2#23, sin 2 2#i3, and the neutrino mass squared 
ratio R y versus the phase parameter old for &d = 8.7. The horizontal lines denote observed 
values sin 2 29f b 2 s = 0.875 ± 0.024, sin 2 29f bs x 10 = 0.98 ± 0.13 and R obs x 10 = 0.3241°;°^. Our 
predicted value for sin 2 2#23 is well satisfied the obtained experimental bound sin 2 2823 s > 0-95. 

by using the input value [11] Am 2 2 = 0.00243 eV 2 . We also predict the effective Majorana 
neutrino mass [12j (m) in the neutrinoless double beta decay as 



(m) = \miU^ + m 2 C/ 2 2 + m 3 J7 e 2 3 | = 0.0034 eV. 



(3.17) 



This value is a magnitude which may be observed by a future neutrinoless double beta decay 
experiment. 

Finally, we list the predicted values of the CKM mixing parameters and down-quark mass 
ratios, although they are essentially the same as those in the previous model [8]: 

\Vus\ = 0.2271, \V cb \ = 0.0394, \V ub \ = 0.00347, \V td \ = 0.00780, (3.18) 

8 q CP = 59.6° (J q = 2.6 x 10~ 5 ), (3.19) 
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m b 



(3.20) 
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0.0569, 
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m d 
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(3.21) 



Here, we have used a d = 25, £0 = 0.0115, and (0i,0 2 ) = (177.0°, 197.4°). The observed values 
are as follows: \V US \ = 0.2252 ± 0.0009, \V cb \ = 0.0409 ± 0.0011, \V ub \ = 0.00415 ± 0.00049, 
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\Vtd\ = 0.0084 ± 0.0006, J = (2.96^) x 10" 5 JO], and rft = 0.045±g;gig ) r£ 3 = 0.060 ± 0.005 



'12 — U - UO,3 -0.003' '23 



0.019 ±0.006 113 . 



4 Concluding remarks 

In conclusion, we have proposed a new neutrino mass matrix form within the framework 
of the yukawaon model, in which we have only two adjustable parameters, (ap, ap), for PMNS 
mixing and neutrino mass ratios. We obtain reasonable results for PMNS mixing and neutrino 
mass ratios as shown in Eqs.(3.13) - (3.17) for the parameter values (od,Q!£)) = (8.7,12°). 
As seen in Fig. 2, it is worthwhile noticing that only when we choose a reasonable value of 
R v ~ 0.033, we obtain a reasonable value of sin 2 2#i3 ~ 0.09. Also, we would like to emphasize 
that our prediction gives a sizable value of (m) — 0.0034 eV in spite of the normal mass hierarchy 
model (in spite of m v \ ~ 0.0006 eV, m u2 ^ 0.00899 eV and m u3 ~ 0.05011 eV). 

Of course, we have also obtained reasonable results for CKM mixing and quark mass ratios 
as same as those in the previous paper [8] . 

The present model gives successful results from the phenomenological point of view. How- 
ever, we still have some of theoretical problems. One of the major problems is why only Yd 
takes the mass matrix form with X 2 (not A3). In Ref.[6j, the form A3 has been understood by a 
symmetry breakdown U(3)xU(3)' — > U(3)xS3. However, for the form A2, the model is still in a 
phenomenological level. We have been able to remove the unnatural term \£ v term in Eq.(1.2)], 
but we still have £q term in the sector. These problems are our future tasks. 

Appendix 

In this Appendix, we discuss a derivation of the mass matrix relations (2.2)-(2.9) from 
superpotential. We assume the following superpotential W = W e + W u + Wr + Wd + W u + Wd'- 
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The VEV matrix relations (2.2) - (2.9) are obtained from SUSY vacuum conditions, dW/dQ^ = 
(A = e,u, ■ ■ ■). Since we assume that all B fields take (G) = 0, SUSY vacuum conditions with 
respect to another fields do not lead to meaningful relations, because such conditions always 
contain, at least, one (0). 

In Eqs.(A.5) and (A. 6), we have introduced fields E", E' d , E u and E d in addition to E" and 
E in order to distinguish the R charges of 6 M and @ d from e . All VEV matrices (E) are given 
by the forms (E) oc 1 as seen in (A. 10). The VEV matrix relations (2.2) - (2.9) have already 
been presented by replacing (E) —> 1. 

We list the assignments of SU(2)ixSU(3) c xU(3)xU(3)' and R charges for additional fields 
in Table 3. The assignments of R charges are done so that the total R charge of the superpotential 
term is R(W) = 2. We have 17 constraints on the R charges of the fields from Eqs.(2.1) and 
(A.l) - (A. 6), while we have 34 fields even except for B fields in Tables 1 and 3. Therefore, 
we cannot uniquely fix R charge assignments of those fields. Here, let us give only typical 
constraints: 

2r x = r" E -f E = r E -f" E = r" Eu - f Eu = r" Ed - r Ed , (A.7) 

2r = f Ye ~r E = r Y D - r E = r Yu + 2r Eu - r Eu = r Y d + 2r Pd - r" Ed . (A8) 

From Eq.(A.7), we obtain r"+f E = r E +f E . When we take R(E")+R(E") = R(E)+R(E) = 
R(P d ) + R{P d ) = 1, we can introduce the following superpotential: 

W e ,p = ^Tr[EEP d P d ] + ^Tr[EE]Tr[P d P d ], (A.9) 

from which, we obtain relations (E)(E) oc 1 and (P d )(P d ) cx 1. We assume following specific 
solutions of those relations: 

±{E) = ±(E} = 1, (A.10) 
VE ve 

— (Pd^ = ^(P d ) = diag(e-^,e-^,l), (All) 
Vp vp 

as the explicit forms of (E), (E) and {P d )- We assume similar superpotential forms for (E,E), 
(E U ,E U ), (E d ,E d ), (E",E"), (E"»,E%), (E'',^) and {E>,E>). 

The term fi d E d in Eq.(A.6) has been introduced in order to adjust the down-quark mass 
ratio m d /m s as seen in Sec. 3. Similar bias terms E and E u cannot appear in the lepton and 
up-quark sectors, because we take R(E) / R(E d ) and R(E U ) ^ (E d ). 
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Table 3: Assignments of SU(2) L xSU(3) c xU(3)xU(3)' and R charges 
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